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Abstract 

Graphite is an example of a layered material that can be bent to form fullerenes which promise 
important applications in electronic nanodevices. The spheroidal geometry of a slightly elliptically 
deformed sphere was used as a possible approach to fullerenes. We assumed that for a small 
deformation the eccentricity of the spheroid e <C 1. We are interested in the elliptically deformed 
fullerenes C70 as well as in Ceo and its spherical generalizations like big C240 and C540 molecules. 
The low- lying electronic levels are described by the Dirac equation in (2+1) dimensions. We show 
how a small deformation of spherical geometry evokes a shift of the electronic spectra compared to 
the sphere. The flux of a monopole field was included inside the surface to describe the fullerenes. 
Both the electronic spectrum of spherical and the shift of spheroidal fullerenes were derived. 
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I. INTRODUCTION 

Fullerene molecules jl| are carbon cages which appear in the vaporization of graphite. 
One of their most beautiful features from a formal point of view is their geometric character 
and the exciting possibility of producing them in all sorts of geometric shapes having as 
building blocks sections of the honeycomb graphite lattice. The most abundant of them is 
the most spherical Cqq molecule. The shape of the Cqq molecule is that of a soccer ball, 
consisting of 12 pentagons and 20 hexagons. However, some fullerenes as C70 are slightly 
elliptically deformed with the shape being more similar to an American football. Fullerenes 
belong to a sort of carbon nanoparticles. 

Carbon nanoparticles, which are expected to have important implications for the devel- 
opment of electronic devices, flat panel displays, nano- switches, etc., have recently received 
great attention of both experimentalists and theorists (see Q]). High flexibility of carbon 
allows one to produce variously shaped carbon nanoparticles: nanotubes, fullerenes, cones, 
toroids, graphitic onions and nanohorns. Particular attention was given to peculiar elec- 
tronic states due to topological defects which were observed in different kinds of carbon 
nanoparticles by scanning tunneling microscopy (STM). For example, STM images with 
five-fold symmetry (due to pentagons in the hexagonal graphitic network) were obtained 
in the Ceo fullerene molecule j^]. The peculiar electronic properties at the ends of carbon 
nanotubes (which include several pentagons) were probed experimentally in 0, sj]. 

By its nature, the pentagon in a graphite sheet is a topological defect. Actually, as was 
mentioned in Ref. (|, fivefold coordinated particles are orientational disclination defects in 
the otherwise sixfold coordinated triangular lattice. The local density of states was found in 
the vicinity of a pentagonal defect for spherical fullerenes [?], s| . Moreover, disclinations are 
generic defects in closed carbon structures, fullerenes and nanotubes, because, in accordance 
with Euler's theorem, these microcrystals can only be formed by having a total disclination 
of An. According to the geometry of the hexagonal network, this implies the presence of 
twelve pentagons (60° disclinations) on the closed hexatic surface. 

Investigation of the electronic structure requires formulating a theoretical model describ- 
ing electrons on arbitrary curved surfaces with disclinations taken into account. An impor- 
tant ingredient of this model can be provided by the self-consistent effective-mass theory 
describing the electron dynamics in the vicinity of an impurity in graphite intercalation 
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compounds The most important fact found in [s| is that the electronic spectrum of a 
single graphite plane linearized around the corners of the hexagonal Brillouin zone coincides 
with that of the Dirac equation in (2+1) dimensions. This finding stimulated a formulation 
of some field-theory models for Dirac fermions on hexatic surfaces to describe electronic 
structure of variously shaped carbon materials: fullerenes Q, [ll| and nanotubes 12]. 

The Dirac equation for massless fermions in three-dimensional space-time in the presence 
of the magnetic field was found to yield N — 1 zero modes in the N-vortex background 
field [l3]. As was shown in Ref. [ijj], the problem of the local electronic structure of fullerene 
is closely related to Jackiw's analysi s |13f] . Notice that the field-theory models for Dirac 
fermions on a plane and on a sphere [15| were invoked to describe variously shaped carbon 
materials. Recently, the importance of the fermion zero modes was discussed in the context 
of high-temperature chiral superconductors and fullerene molecules. 

The most spherical fullerene is the Cqq molecule nicknamed a 'bucky ball'. Others are 
either slightly (as C70 whose shape is more like an elliptic deformation) or remarkably de- 
formed. We are interested here in the Ceo molecule as well as in its spherical generalizations 
like big C240 and C540 molecules with the symmetry group of the icosahedron, and also in 
the elliptically deformed fullerene C70 and its relatives. Big fullerenes are used to store 



radioactive material and inhibit enzymes related to different viruses 16l. Il7l|. 



II. THE MODEL 

Almost all fullerenes are only slightly elliptically deformed spherical molecules, e.g., C70 
and its relatives. We start with introducing spheroidal coordinates and writing down the 
Dirac operator for free massless fermions on the Riemann spheroid S 2 . Pi-molecular orbitals 
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erenes as a free electron model (electron gas) bound on the surface of a sphere were used 
] . We generalize that work to obtain an electronic spectrum for spherical and spheroidal 
geometries with and without the monopole field. The peculiarities of the electronic spectra 
for these two slightly different types of geometries are shown. 

To incorporate fermions on the curved background, we need a set of orthonormal frames 
{e a }, which yield the same metric, g^, related to each other by the local SO(2) rotation, 

e a ^e' a = A%, AieSO{2). 
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It then follows that g^ u = e°e^5 a p where is the zweibein, with the orthonormal frame 
indices being a, (3 — {1,2}, and the coordinate indices /i, v = {1,2}. As usual, to ensure 
that physical observables are independent of a particular choice of the zweibein fields, a local 
so(2) valued gauge field is to be introduced. The gauge field of the local Lorentz group 
is known as a spin connection. For a theory to be self-consistent, the zweibein fields must 
be chosen to be covariantly constant [loj] 

V,e a v := «V« - J%el + (u^e? = 0, 

which determines the spin connection coefficients explicitly 

(a*,)* = e" V*. (1) 
Finally, the Dirac equation on a surface £ in the presence of the magnetic monopole field 

i-fe^^-iA^ = E^ (2) 



is written as 



where V M = <9 M + f2 M with 



n /l = iw a /[7«,7/»], (3) 



being the spin connection term in the spinor representation. 
The elliptically deformed sphere or a spheroid 



l2 +M-L (4) 



a 2 a 2 c 2 



may be parameterized by two spherical angles q 1 = 9, q 2 = <fi that are related to the 
Cartesian coordinates x, y ', z cis follows 

x = a sin 9 cos 0; y = a sin # sin 0; z = c cos#. (5) 

We have assumed that the eccentricity of the spheroid is e <C 1 which in the case c < a 
gives expressions e = ^/l — (^) 2 <C 1 . The metric tensor and the natural diagonal zweibein 
different from zero on the spheroid are 

Q<$>4> = a 2 sin 2 9; ggo = a 2 cos 2 9 + c 2 sin 2 9, (6) 

where a, c > 0, < 9 < ir, < <fi < 2n and 



which, in view of Eq. ([I]) gives the spin connection coefficients 

a 



^ = -u 2 } 



--: 2u. 



(8) 



y/a 2 + c 2 tan 2 6 

In 2D the Dirac matrices can be chosen to be the Pauli matrices, 7 1 = — a 2 , 7 2 = a 1 ; Eq. (J3]) 
then reduces to 

VL^ = iua 3 . (9) 

We have assumed that Ag = and only the monopole field is different from zero. 

The eigenfunctions of the Dirac operator are two-component spinors that satisfy the 
eigenvalue equation 

(10) 



/T 'ax{B,<j>)\ =x ( a\{0,4>) 



A(M)/ \ft(M). 

This system of first order partial differential equations in a and f3 allows separation of 
variables therefore we can isolate the ^-dependence by expanding the spinors into the Fourier 
series 

if, fa ^ f a x(°^)\ ^ exp im(j) f a Xm (6)\ 1 3 , , 

where m are half-integers since we work with the spin ~ field. Then the general form of the 
Dirac equation, Eq. (j2j), on the spheroid becomes 



d e + \ I cot 2 9 + ( - a Ym + (l - 2A+ J 1 + tan 2 e \ ^ 



a 2 cos 2 6 + c 2 sin 2 6 Ea Xm {6), 



do 



< C ot 2 e +{ c)m + ll + 2AJl + (l) tan 2 ,)^ 



x/a 2 cos 2 6 + c 2 sin 2 9 E(3 Xm (9). (12) 

The number m may be called the projection of angular momentum onto the polar axis. If 
a = c = R, where R is the radius of a sphere, Eq. (fT2l) becomes the Dirac equation for 
sphere geometry. 

Now we want to find an electronic spectrum for the sphere and spheroid analytically and 
numerically, respectively; therefore, we firstly assume that pentagon defects represented in 
this model by the monopole field A^ = 0. Next, we want to separate the equations for the 
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spinor components a and (3. This can be done by taking the square A (Laplace operator) 
of the Dirac operator V for spheroidal geometry. Finally, we find the equations 



n d e sm9d e +\ cot 2 6+ 
sin 9 \ \ </ 



m 2 -a 3 



cot9 



rn 



'cot 2 fl+ 



2 sin 9 4 sin 9 



a\ m (9) 



{a 2 cos 2 9 + c 2 sin 2 9) E 2 



OC\ m {9) 



(13) 



Further simplifications come from the change of variables x = cos 9, x G [—1,1], which 
converts Eq. fflBl into the generalized hypergeometric equations 



d 

dx 



i — 

dx 



x 



or 



m 2 +a 3 



x 



1 



x 2 4(1 - x 2 ) 



^ Am \pC) 



a 2 x 2 + c 2 (l -x 2 )) E - 1/4 



(14) 



The replacement x — > —x (or m — > —m) is equivalent to changing a for f3. Thus, 
the upper and lower spinor components are conjugate with respect to mirror reflection. 
Equation (ff4l) is singular at the poles of the spheroid x — ±1. We redefine the unknowns 



(1 _ x )V2|m-l/2| (1 + a; )V2|m+l/2|^ m ( :!; ) 
(1 - x)l/2[r»*+V2|(l + Xy/^-WrjxJix) 



(15) 



and use that in our model of a slightly deformed sphere c ~ a; hence - = 1 ± 8, where 5^1 
is small deformation of the sphere. Neglecting the second and higher order powers of S we 
solve Eq. ffT4l) to the first order in 5 and by using redefinition ffT5|) we arrive at the separate 
equations of hypergeometric type in £ Am and rj\ m 



2(m — l)x) 



d 

dx 



m{m - 1) + A 2 - 1/4 + f(x) }t Xm (x) = 0, (16) 



where 



f(x) = 5(T2m 2 ± 2a 2 E 2 =f mx =f 2a 2 E 2 x 2 ); X = Ea 



For the above calculations we have assumed that m — 1/2 < 0. A similar solution can be 
also found for the case m — 1/2 > 0. Following the calculations above we can get also 



the equation for function r]\ m (x). T 
spheroid from that for a sphere, see 



re function f(x) is the deviation of the solution for a 



21], and can be perceived as next energy part, i.e. the 
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energy shift for a spheroid compared to a sphere geometry. Thus the expression f(x) can 
also be called as a perturbation part of Equations f|T6|) . 

To find the spectrum for spherical geometry a = c = R, we have to put the expression 
f(x) = 0, so the case if 5 = 0. The spectrum can be found in the form 



Sphere ={n+\m\ + l/2)\ 



(17) 



with non-negative integer n > 0, with n being the order of Jacobi polynomials, see 2l|]. The 
eigenvalues A for the sphere of the unit radius S 2 are nonzero integers 



A 



sphere 



±1,±2. 



(18) 



and indeed the Dirac operator has no zero-modes. The spectra of the spherical geometry as 
the numerical calculations of Eq. (|T2|) with = and for a = c are illustrated in Fig. 1 
and fit the analytical results in (18). 




FIG. 1: The electronic spectra of spherical geometry |A| sp /j ere , where a = c. 

To find the electronic spectrum also in the case of spheroidal geometry when f(x) ^ 0, 
we solve Eq. (TT2]) for two cases, when a < c and a > c, i.e. for two different types of an 
elliptically deformed sphere. The numerical results are shown in Fig. 2. 

As can be seen from Fig. 2, the spectra display the energetic shifts compared to the 
spectra of the sphere in Fig. 1. The shift increases or reduces the initial electronic spectrum 
of sphere depending on the type of elliptic deformation. The shift is bigger with increasing 
value of the modes of the electronic spectra, which can be also seen from the structure of 
f(x) in Eq. ( JT6l) . The perturbation parameter 5 = 0.3 was used to obtain the electronic 
spectra of the Dirac equation for the spheroidal geometry. The angular momentum m — 1/2 
was used in the calculations. All the spectra are mirror symmetric with respect to the y-axis. 
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FIG. 2: The electronic spectra of spheroidal geometry \\\ sp heroid where, a < c and a > c (gohi£ 
from left to right). 

III. DIRAC EQUATIONS FOR SPHEROID WITH MONOPOLE FIELD 



We assume that the parameter of perturbation and the eccentricity of the spheroid e <C 1, 
so we can use the magnetic monopole field inside the surface to obtain CVo fullerene or bigger 
fullerene molecules like C240 and C540 also with a small elliptical deformation. The area of 
surface for a small elliptically deformed sphere, spheroid, e.g. for oblate spheroid (a = b) 



can be formulated 



22] 



S « 2n(a 2 + c 



,arctanh(e) 



)■ 



In the case of small eccentricity the surface area of a spheroid becomes the surface area 
of a sphere S ~ 2-ko 2 . It means that we can include the existence of a fictitious magnetic 
monopole charge g inside the surface of the spheroid with the structure as in Ref. 15]. 
The values of g, e.g. for tetrahedron and icosahedron structures required 1/2 and 3/2, 
respectively. With the monopole field taken as Ag = and A^ = ^cos9, where j/2 = g 
(for structures above j = 1,3), the resulting Dirac equation (fl2l) for the spheroid with the 
monopole field reads 



d e + \ cot 2 9 



m 



+ ( 1 - 3 -\la 2 cos 2 9 + c 2 sin 2 9 



cot 9 



P\m(9) 



a 2 cos 2 9 + c 2 sm 2 9 Ea Xm (9), 



do 



'cot 2 9 + ( - 1 /;/ 
a 



J 



1 + -Va 2 cos 2 fl + c 2 sin 2 # , 

a 2 



cot 9 



a\ m {9) 



Va 2 cos 2 9 + c 2 sin 2 9 E(3 Xm (9). 



(19) 



Unfortunately, general solutions to Eq. (|T9j) are not available yet for spherical and spheroidal 
fullerenes, so we do not have initial conditions for numerical calculations which are very 
sensitive to them. However, we will present some analytical predictions of the electronic 
spectra for spherical and spheroidal fullerenes from the square of Dirac operator in Eq. ffT^l) . 

To find an analytical expression for the shift of the spheroidal fullerenes, we put for 
simplicity in Eq. (|T9j) the value of the angular momentum m = 0. For elliptically deformed 
fullerenes in the special case where m = with using the model of a slightly deformed sphere 
where - = 1 ± 8, the square of Dirac operator (19) was found in the form 



1 ^u,^ 5 



4 sin 2 9 
(A 2 + (j/2) 2 )(l±25sin 2 



a 3 ) 2 , ±j6 cos 2 6 
— — + a 3 



V(l ±2<5sii 



a\ m (0) 



1/4 



(20) 



sin 9 i mii in % / 1 i -+- \o sirr 

where 5 <C 1 is small deformation of the sphere. We neglect the second and higher order 
powers of S in the calculations. 

If the perturbation parameter 5 = 0, we get the square of Dirac operator for spherical 



fullerenes with the electronic spectrum as in Ref. [Tl|. To get connection with Section II, 
the spectrum with the monopole field for the spherical fullerenes can be rewrite from the 
right hand side of Eq. (120]) to the form 

(|A| sphere -l/2) 2 -/, (21) 



X 2 

A field 



with |A| S pfe ere — 1/2 = n + \m\, see [21]. The value of the electronic spectrum X field is shifted 



(decreases) by the value of the charge g of the monopole field. Moreover, the presence of the 
monopole restricts possible values of the angular momentum [is| . so that m >\\ j \ — 1 \ /2 
and, therefore, the values of m change contrary to structures without a monopole field inside, 
as in Section II. It means that for fullerenes the angular momentum m can get the values 
m = 0, 1.... The spectra of fullerenes are appended by the monopole charge compared to 
Eq. f|T3|) for spherical geometry. So the spectra of fullerenes are greatly dependent on the 
value of the monopole field. In the case of tetrahedron and icosahedron structures expression 
(21) can predict the existence of zero modes where X field = 0. 

If S ^ 0, the shift of the spheroidal fullerenes from the spherical ones in Eq. f)20p were 
found. Moreover, when we use the substitution x = cos 9 as in Section II, we can rewrite 
the shift on the right hand side of Eq. (120]) in the following form 

f(x) = 25(X 2 + g 2 )(±lTx 2 ); X = aE. (22) 



The electronic spectrum will be shifted by the value ±2<5(A 2 + g 2 ). The function f(x) is 
the deviation of spheroidal from spherical fullerenes, similarly to the previous section, and 
can also be perceived as the energy shift to spheroidal fullerenes compared to spherical ones. 
Moreover, the perturbation function f(x) is in the case when m = and g = the same as 
the shift in Eq. ( I16p for spheroidal geometry. So we can expect a similar behavior (shifts) 
of the spectra of spheroidal fullerenes (decreases by the value of the monopole charge g) as 
in the case of spheroidal geometry without a monopole field, see Fig. 2. 



IV. CONCLUSION 



To find the electronic spectra of the CVo fullerene and its relatives, we have used the 
model of a slightly elliptically deformed sphere, the spheroidal geometry (5^0), as distinct 
from the Cqq fullerene where the spherical geometry (5 = 0) was used. The Dirac equation 
in (2+1) dimensions for slightly elliptically deformed fullerenes with monopole field inside 
the surface was evaluated. The discrete spectrum of energy for both types of geometries 
was found. The electronic spectrum and the shift of the spheroidal geometry in Eq. (|T6|) 
contrary to the sphere was calculated both analytically and numerically. Figure 2 shows the 
shift of the spectra for spheroidal geometry. In the case of spherical and spheroidal fullerenes 
the electronic spectrum and shift were derived analytically. The spectra of spherical and 
spheroidal fullerenes decreasing by the value of monopole charge g were found. The expres- 
sion / (x) for deviation of the solution for a spheroid from that for a sphere in Eq. (|T6|) was 
found the same as for the deviation between spheroidal and spherical fullerenes in expression 
(22). So we can expect a similar shift of the spectra for real elliptically deformed fullerenes 
when the magnetic monopole field has to be included inside the surface to simulate pentagon 
defects and create fullerenes. Zero energy modes dictated by the charge g were found for 
fullerenes contrary to spherical geometry without the monopole field inside. The shift of the 
electronic spectra of spheroidal, in contrast to spherical fullerenes, gives rise to reduction or 
increase in the conduction bandwidth depending on the type of elliptical deformation (a < c, 
a > c) . Due to this, the crystals made of these deformed fullerene molecules, when doped 
should be poorer or better conductors than the spherical ones. Moreover, the spherical 
fullerenes as Cqq are stable towards fragmentation than the other bigger fullerenes [23[ , fol- 
lowing our analysis, also with small elliptical deformations and with a shift of the electronic 
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spectra. The presence of the magnetic monopole field with the charge leads to a decrease in 
the electronic spectra and shifts for spherical and spheroidal molecules, respectively. The de- 
crease is smaller with increasing value of the modes of the electronic spectra and, therefore, 
for low-lying electronic levels the spectra of the spheroidal fullerenes Cjq or C240 and C540 
could be shifted to a lower magnitude. The very big structures like C 960 and Ci 50 o become 
more deformed, faceted and can no longer form a free-electron model like the electronic 



shell 18| . which was the assumption for this model. For these structures the deviation from 
the sphericity is larger when the pentagon defects are localized at the opposite poles. In 
the case when the poles are far away from each other we obtain the structure of nanotubes, 
and for the exact description some new model related to that proposed here should be used. 
We hope that the knowledge of the shifts of the electronic spectra of spheroids could be 
useful for experimentalists for choosing the optimal energetic scale for different types of 
fullerenes. Finally, we think that the spheroidal geometry approach could also be related to 
other physical problems with slightly deformed spherical structures that are common in the 
nature. 
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